The effect of perception anisotropy on particle systems 
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We consider a microscopic model (a system of self-propelled particles) to study the behaviour of 
a large group of pedestrians walking in a corridor. Our point of interest is the effect of anisotropic 
interactions on the global behaviour of the crowd. The anisotropy we have in mind reflects the 
fact that people do not perceive (i.e. see, hear, feel or smell) their environment equally well in 
all directions. The dynamics of the individuals in our model follow from a system of Newton-like 
equations in the overdamped limit. The instantaneous velocity is modelled in such a way that it 
accounts for the angle under which an individual perceives another individual. 

We investigate the effects of this perception anisotropy by means of extensive simulations, very much 
in the spirit of molecular dynamics. We define a number of characteristic quantifiers (including the 
polarization index and Morisita index) that serve as measures for e.g. organization and clustering, 
and we use these to investigate the influence of anisotropy on the global behaviour of the crowd. 
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I. INTRODUCTION 

During the last two, three decades, the field of crowd 
dynamics has emerged as the natural sciences' reaction 
to questions arising from social sciences, population 
biology and urban planning. The roots and philosophy 
of crowd dynamics are very much in the spirit of 
statistical mechanics, molecular dynamics, interacting 
particle systems methods and the theory of granular 
matter, as such treating individual humans nearly as 
non-living material (cf. e.g. [H 03 QUI HU [H] or the nice 
overview [Fj and references cited therein). A justification 
for this approach lies in the fact that the individuals' 
personal will is more or less averaged out if one looks at 
the crowd as a whole. From this perspective, it can be 
considered as (stochastic) noise, superimposed on some 
'clean' (deterministic) dynamics. 
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An evident and important difference between peo- 
ple and molecules or grains is the fact that people clearly 
have own opinions, own irritations, front and back sides, 
etc. Our degree of perceiving our surroundings highly 
depends on the direction of looking. We mainly base 
our walking behaviour on what we see, and clearly what 
happens in front of us thus has more influence than 
what happens behind us. A modification or extension 
of physics-inspired models is needed to incorporate this 
kind of anisotropy in the interactions between individu- 
als. This paper investigates the effect of anisotropy on 
the global behaviour of a group of pedestrians. 

Our focus is on the simulation of a scenario where 
pedestrians move in a long corridor. We might relate 
this situation to evacuation of people from a building 
(cf. [HHi]). I* * s a sane to assume tna t, these evacuees 
have an intrinsic drive to move towards the exit (i.e. 
one side of the corridor), and moreover that there view 
is focused in the same direction. Investigating the effect 
of anisotropy on the large-scale behaviour of the crowd 
therefore relates to assessing the escape process. 
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In Section |TT] of this paper the model is presented 
and explained. Section [III| is the main part of the paper. 
It describes the exact scenario of our simulations and 
the definitions of the quantities we use for assessing the 
results (polarization index, projected density, Morisita 
index). Moreover, in this section the simulation results 
are presented and discussed. Conclusions and an outlook 



on possible future work are given in Section IV 



II. A MODEL FOR ANISOTROPIC 
INTERACTIONS BETWEEN PEDESTRIANS 

We represent pedestrians by point particles having 
masses m; . They are located in a long corridor of length 
L and width B. Here, the word 'long' refers to the fact 
that at the time scales we focus on, the pedestrians are 
not able to reach the end of the corridor. Interactions 
between pedestrians are short-ranged. We therefore sup- 
pose that the correlation length in the system is less than 
or equal to a certain L <C L and we can subdivide the 
corridor in an array of rectangles (width B and length L) , 
which are all duplicates of each other. We thus have a 
scenario with periodic boundary conditions. Our domain 
of interest is therefore a rectangular box 
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2' r 
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(1) 



with periodic boundary conditions in one direction and 
impermeable walls in the other direction. The corridor 
contains N < oo pedestrians. For all i E {1, ...,N} and 
t > 0, the vector fi(t) = (xi(t) , Ui{t)) € CI represents the 
position of the i-th pedestrian at time t. We denote its 
velocity by Vi(t). 

We assume that the governing equation of motion 



Tdrive 



-(Vi(t)-v des ) = F l soc + F 1 



phys 



(2) 



The equation describes the motion of the i-th individual, 
which has mass m; and which moves with velocity Vi(t). 
However, he/she tries to move according to its desired 
velocity Vdes- Here, Tdrive is the characteristic relaxation 
time related to attaining the desired velocity. Its actual 
velocity is moreover perturbed by two 'forces'. The 
word 'force' is used since Q can be regarded as an 
overdamped limit of a Newton-like equation (cf. [1] for 
this Newton-like way of modelling). There is a physical 
force p"P hys that acts on the individual to describe the 
effect of the non-living environment (geometry) . In this 
paper we only take into account the influence of walls 



on pedestrians, that is: F- 



phys 



Fj 



Furthermore, 



pedestrian i experiences a so-called social force F? oc due 
to the presence of other individuals, which influences the 
motion of this particular pedestrian i. 



Individuals are influenced by the walls as soon as 
they come too close, say within a distance R wa n- We 
model these impermeable walls by means of a strong 
repulsive force F™ al1 acting on pedestrian i: 



ffwall TP f I Rwall 



(3) 



where n is the unit normal pointing from the corre- 
sponding wall into the corridor, Fwaii is the strength 
of the repulsive force and d is the distance to the wall 
for pedestrian i. The word 'strong' here implies that 
this force is not just a contact force, but has a longer 
range. Typically, this makes individuals avoid walls 
before touching them. 

Furthermore, very much in the spirit of pQ, we 
specify the social force by 



(4) 



where: 



• Cli is the collection of the position vectors of all in- 
dividuals which are within a distance R cu t to pedes- 
trian i. In other words, pedestrians interact only 
when they are close enough to each other; 

• we assume that the interaction potential W de- 
pends only on the relative position of the two pedes- 
trians i and j and not on their relative velocity. 




FIG. 1. Schematic drawing of the variables ri, fj, fj — f, 
involved angle 6 riTj and unit vectors e x , e y . The vector Vd es 
can in principle have arbitrary direction. For the sake of clar- 
ity in the picture, it is chosen to be parallel to e x . Here, 6r irj 
is the angle under which an individual positioned in fi and 
heading along Vdes, perceives location fj. 



Specifically, W takes the form, 

W(n - r,) = U{\n - f^])(|(l + <y) + \{l- <t)cos 9 rtTj ). 

(5) 

The parameter a € [0, 1] is called potential of anisotropy, 
for which a — means that the anisotropy effects are 
maximal, and a — 1 means that the potential is isotropic. 
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The angle of perception 9 rir . is the angle under which an 
individual positioned in fi, while moving in the direction 
of Vdes, perceives location fj, see Fig. [T] The precise type 
of interaction is hidden in the structure of the function 
U . In particular, we distinguish between two types of 
interactions, namely 

• only repulsive interaction (for simplicity denoted by 
R case ). Pedestrians repel each other when their 
separation distance is smaller than R R (called ra- 
dius of repulsion in R case ) and do not interact at 
larger distances; 

• both repulsive and attractive interaction (denoted 
by AR case ). For this type of interaction, indi- 
viduals repel each other when their separation dis- 
tance is smaller than R^ R (called repulsive radius 
in AR case ). However, when they are separated by 
a distance between R^ R and R£ r (called attractive 
radius in AR case), they are attracted to one an- 
other. They do not interact outside these regions. 

For these two cases R cut = R R , respectively R cut — Ra R - 
The second case is an extension of the first one. The 
R case we might regard as a population of individualis- 
tic people that simply try to avoid each other. In the 
AR case there is also some social cohesion, as they try 
to keep the group together. In the following, we describe 
the precise structure of our potentials. 



2. AR case 

If we want to involve a both attractive and repulsive 
way of interactions, then we take 



U(s) 



where 



U(s) 



F AR (t 

U(s) 




R AR ln 



R£ R ' 



C 2 tfs<Rf R ; 

if R AR <s<R AR ; 
if a > R* R ; 

(7) 



F AR (j - (R? R + R AR )s 2 + R AR R AR s) 
R*R{RA R - RAR) 



Ci- 

(8) 

In the above, the constants C\ and Ci are such that U 
is continuous. Fig. [3] shows an example of this kind of 
interaction potentials. 

U(s) 




FIG. 3. Typical example of the interaction potential as given 
in @-©- 



1. R case 

In the repulsive case we take 

{ F R (s - R R - R R In — ) if s < R R - 
U(s) = \ { r r R R> r ' (6) 

to if s > R R . 

See Fig. [2] for an example of an interaction potential of 
the above form. 



U(s) 




s 



FIG. 2. Typical example of the interaction potential as given 
in J6j. 



III. SIMULATION: SET-UP AND RESULTS 

For all simulations, the number of pedestrians N is 
an integer multiple of 10. Initially, we place these in- 
dividuals on a lattice of N x = ^ rows of 10 pedestri- 
ans each. The distance between two succeeding rows is 
always Ax = jr. Note that there is an N here, not 
N x : the pedestrians are thus distributed initially in a 
domain about one tenth the length of the total corri- 
dor (see Fig. [4]). The distance between two pedestri- 
ans in the same row is Ay — A. Let (h,yj) denote 
the position of pedestrian k := i + 10(j — 1), where 
i = 1, . . . , 10 and j — 1, . . . , N x . Then the coordinates 
are k = =^ + Ax{i - §) and y 3 , = =f + Ay(j - \). We 
use these initial conditions for all simulations. 

The simulation time corresponds with a real life time 
Tobs — 100 s. This time is large enough to witness a 
stable profile for all simulations. We give in Table [I] a 
summary of the parameters and their values used in both 
the 'repulsive' and 'attractive and repulsive' interaction 
potentials. 

For more technical details about this kind of simula- 
tions, the reader is referred e.g. to [THUS]. 
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FIG. 4. A snapshot of the positions of the 100 pedestrians at 
time t = s. This initial configuration is used in all simu- 
lations. The markers that indicate the individuals' positions 
are smaller than in Figs. [5^10| just to avoid overlap here. 




FIG. 5. A snapshot of the positions of the 100 pedestrians at 
time t = 100 s, plot for the R potential with a = 0.0 (most 
anisotropic interactions) . 



To illustrate simulation results, in Figs. [5|[10| snap- 
shots of the position of N — 100 pedestrians at time 
t — 100 s are shown in both the R case and the AR case 
for a = 1.0, er = 0.5 and a = 1.0 (fully isotropic), 
respectively. 

Excepting Fig. [8j all figures depict a completely 
symmetric profile with respect to the line y = 0. This is 
a natural result of the symmetry in the initial data and 
symmetry w.r.t. the direction of Vdes in the interactions 
[This is just a property of the hyperbolic conservation 
laws to govern the mean- field behaviour [IS].]. Note 
also the strong tendency of the system to maintain 



TABLE I. Model parameters, symbols, units and reference 
values. We do simulations in two cases. The AR case : pedes- 
trians interact in a both attractive and repulsive way and the 
R case : pedestrians only interact in a repulsive way. NB: M 
and S x ■ S y relate to the Morisita index and are specified in 
Section lDTA~3l 
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Unit 
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20, 40, 60, 80, 100 
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FIG. 6. A snapshot of the positions of the 100 pedestrians at 
time t — 100 s, plot for the R potential with a = 0.5. 
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FIG. 7. A snapshot of the positions of the 100 pedestrians at 
time t = 100 s, plot for the R potential with a — 1.0 (isotropic 
interactions) . 
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FIG. 8. A snapshot of the positions of the 100 pedestrians at 
time t = 100 s, plot for the AR potential with a = 0.0 (most 
anisotropic interactions). Note that the symmetry is broken 
due to discretization errors. 
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FIG. 9. A snapshot of the positions of the 100 pedestrians at 
time t — 100 s, plot for the AR potential with a = 0.5. 




x (m) 

FIG. 10. A snapshot of the positions of the 100 pedestrians 
at time t = 100 s, plot for the AR potential with a — 1.0 
(isotropic interactions). See main text for the used initial 
conditions. 



(or produce) organized patterns. These issues will be 
addressed in more detail in Section IIIIBI 
Fig. [8] is an exception in the sense that the distribution 
of the crowd is not symmetric around the line y = 0. 
This cannot be explained from the model equations, as 
these imply that the distribution should be symmetric. 
This effect must be due to discretization errors and their 
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propagation with respect to time. 

In the AR case for a = 1.0 (cf. Fig. 10 ) we see that the 



people move in a crystal-like formation parallel to the e x 
- axis, which remotely resembles the 2D crystallization 
patterns at low temperature pointed out in [7J. The 
group as a whole is very compact. In the sequel we will 
use the words clusters/clustering for a situation in which 
pedestrians move close together. 

However, if we compare Fig. [10] to Figs. § and [9j we 
observe that introducing anisotropy by setting a < 1 
decreases the amount of structure and clustering. In the 
R case , all three values of a allow for a well-structured 
way of moving, but there is no clear clustering. It is 
evident from the figures that the people are spread over 
the whole corridor. 

Now, we introduce a number of measurable quanti- 
ties to help quantify the above statements. 



where (v(t) ■ e$) n denotes averaging over the total 
number of N individuals: 



7(t) ■ e ( ) 



N 



N 



, Z e {x,y}; 



(13) 



\\6i(t),0(t)\\ := min feez |0;(t) - 6{t) + 2kn\ denotes 
the angle in [0,7r] between Vi(t) and the average 
direction of motion. 



2. Projection of the pedestrian density on the e x axis 

To examine the distribution of our crowd in the di- 
rection that corresponds to the desired velocity, we can 
consider the number of pedestrians that are located in 



S' 



B 



D 



{(x,y) < y < —, \r) - x\ < e} 



A. Definitions of the measured quantities 

1. Polarization index 

Inspired by [B], we define the (time-dependent) polar- 
ization index p of a group of people as the average angu- 
lar deviation from the mean propagation direction. Here, 
the average is taken over all individuals. Note that zero 
polarization means that all the people move in parallel. 
This definition in words translates into 



p(t) :- 



1 N 



(9) 



for the instantaneous polarization. We are also interested 
in the time average of p, which is defined as 



(10) 



In the above we use that 

• (q) t := —^j-^- denotes the time average of an 
arbitrary quantity q, based on the values q(tj) at 
time tj, j G {1,2, M}; 

• 0i(t) G (— 7r,7r] is the direction of motion of the 
pedestrian i at time t. It is defined to be the angle 
0i(t) such that 



tan(0i(t)) 



Vj(t) ■ e y _ 
Vi(t) ■ e x ' 



(11) 



0(t) G (~7i", 7r] is the mean direction of motion of 
the pedestrians group. It is formally defined as the 
angle such that 



tan(0(t)) 



(v(t) ■ e 2 )jy 
(v(t) ■ ei)jv' 



(12) 



Here, < e <C 1, and denotes a narrow strip paral- 
lel to the e y axis, centered at x-position 77 and of width 
i]. We could simply plot the number of people in 
as a function of 77. Since pedestrians are point masses 
in our model, this procedure would produce a discon- 
tinuous, histogram-like graph. Moreover, the number of 
discontinuities would highly depend on the value of e. 
In order to smoothen the results, the individuals' coor- 
dinates have been projected on their corresponding x- 
coordinates. Next, we assign an 'induced density' i to 
each individual i, which can be considered as a mollified 



Dirac delta distribution, see Fig. 11 It was constructed 



in such a way that it has support of width fy. 
The total density in a point 77 € [— 4, |r] is obtained by 




-202 
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FIG. 11. The induced density on the e x axis pg^ . prescribed 
by an individual as a function of the spatial coordinate 77, 
relative to its horizontal coordinate a;,. 

adding all individual contributions: 



A 



: =^2p^,i{v), 



(14) 



6 



where the induced density pg m ,i is given by, 



if \xi - t)\ < ^; 



10 

BL 

geS^i-j)! 2 3BL 11 32 - I" 1 "' '/I — 8' 



otherwise. 



Note that this density distribution is normalized: 



(15) 



N .h N 

pgdn = J] b / 2 pg. ti {x)dx = J2 1 = N - ( 16 ) 



i=l 



i=l 



3. Morisita index 

We subdivide the domain fHn equally sized rectan- 
gular boxes of dimensions S x - S y . Following [B] , we define 
the Morisita index J^vi , which is a measure for the degree 
of clustering in our crowd. The index is the probabil- 
ity of having any two pedestrians together in any of the 
A4 boxes, divided by this probability had the individuals 
been distributed uniformly: 



I M ■= M 



E^i n t (n t - 1) 
N(N — 1) : 



(17) 



where n.j equals the number of pedestrians in box i. See 
Fig. [12] for an example that illustrates the use of the 
Morisita index. The figure also shows that the value of 
the Morisita index depends on the number of boxes M. 
and their distribution. 



B 
2 



FIG. 12. Illustration of the Morisita index. In this example, 
the corridor is subdivided into 32 boxes of dimensions S x -S y = 
^ • X- an d 5 pedestrians are distributed over the boxes. This 
results in a Morisita index of 2£ . 



B. Measured quantities: results 

1. Polarization 

To show the kind of information we can deduce 
from the polarization index, we start by examining the 
evolution of p(t) for N — 100. The combined results are 
shown in Figs. fl3] and[T4) 

In each of the six cases we identify a relatively short 
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FIG. 13. The instantaneous polarization p as a function of 
time. Results in the R case for N = 100, for several different 
values of a. 
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FIG. 14. The instantaneous polarization p as a function of 
time. Results in the AR case for N = 100, for several different 
values of a. 



period of time just after t = 0, during which the polar- 
ization decreases rapidly. Afterwards, there is a state in 
which p fluctuates around a certain average level; this 
state might be an equilibrium. For each a in the R case , 
and for a = 1.0 in the AR case , this average level is 
zero (we needed to zoom in in Fig. 13 to verify that 



these curves really decay below the level as in Fig. 14 1. 
The scenarios a = 0.0 and a — 0.5 in the AR case are 
different in the sense that they do not decay to zero. 
Moreover, if one would zoom in, one would see that the 
oscillations in the R case are less rapid and of smaller 
amplitude than in the AR case . This holds especially for 
a = 0.5 and a = 1.0. Comparing R to AR requires some 
extra care, however. It is difficult to compare them in a 
fair way, because of their intrinsically different nature, 
and because of the (in)compatibility of the tested values 
for the interaction radii (cf. Table |l| . 

Let us focus on the initial rapid decay of p: this 
suggests that the initial configuration is not a favourable 
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state for the system to be in. An immediate relaxation 
takes place, implying spreading of the individuals in 
all possible directions until a more preferable situation 
is reached. Figs. [5] [9]support this statement. Most 
clearly in Figs. [8] and [9] the system evolved away from 
the initial configuration. The same statement is true for 
Figs. [5j [6] and[7J although this is less evident. Once we 
realize however that the particles move in what seem to 
be six horizontal rows, we must indeed conclude that 
the particles have deviated from the initial situation 
in which there were ten rows. Fig. [TU] is somewhat 
different, even though we recognize relaxation, as the 
initial occupation was only about one tenth of the 
corridor length. Spreading in a direction parallel to the 
mean direction of motion does not explain the peak in 
p just after t = however. This is because fluctuations 
in the magnitude of the velocity (that is: the speed) do 
not affect the polarization if all individuals move in the 
same direction. The peak shows that there must have 
been some vertical displacement too. 

We remark that in the AR case p does not neces- 
sarily decay to zero, while it does (or at least: seems 
to do up to fluctuations and noise) in the R case . 
Regarding Figs. [5]-p~0j we can distinguish the graphs 
for a = 0.0 and a = 0.5 in the AR case from the four 
others, since they do not seem to possess the degree 
of order and structure that the other graphs do have. 
Strikingly, these are exactly the cases in Fig. [l4| where p 
does not tend to zero. From this we conclude that there 
is a strong relation between the polarization tending 
to zero, and the preservation or favouring of patterns 
and organization in the system. The cases in which p 
oscillates around a non-zero average in the long run, are 
exactly those in which the initial ordered configuration 
has disappeared after some time. Moreover, the fact 
that the polarization remains positive is an indicator 
that the configurations in Figs. [8] and [9] are not stable. 
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FIG. 15. The time average of the polarization P as a function 
of the number of pedestrians N. Results in the R case , for 
several different values of a. 
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FIG. 16. The time average of the polarization P as a function 
of the number of pedestrians N. Results in the AR case , for 
several different values of a. 



P is shown. For each N, we recognize in P the same 
ordering with respect to a as the ordering we have seen 
before in p (just for N = 100, cf. Figs, p] and p}. The 
ordering between a — 0.5 and a = 1.0 ( R case ) is even 
much clearer when considering the time average. 
Note that, comparing R and AR, the absolute differences 
in P are much smaller than those in p. This is mainly 
because the p-curves in the AR case have smaller range 
but larger domain where they are non-negligible. 



2. Projection of the pedestrian density on the e x axis 



Figs. 17 ( R case ) and [l8| ( AR case ) show the projected 
mollified density that was introduced in Section [ill A 2| 
The graphs show that in the isotropic case (a = 1.0) the 
density profile is completely symmetric around the cen- 
ter of mass. Note that this is only the case if the initial 
conditions are symmetric. For a = 0.0 and a = 0.5 in 
the AR case symmetry is no longer present. This is per- 
fectly sane, since anisotropy (i.e. a < 1) was introduced 
to incorporate asymmetric interactions in our model: a 
pedestrian is more influenced by an other individual in 
front of him than by one behind him. The graphs corre- 
sponding to <7 = 0.0 and a = 0.5 in the R case are not 
given. They possess the same oscillatory behaviour as 
Fig. |17| but without being symmetric around the center 
of mass. They do not provide any further information or 
insight, and thus are omitted. 

The highly oscillatory behaviour in Fig. [17] and 
the omitted graphs reflects the ordered, lattice-like, 
structures we already observed in Figs. [5}{7j Fol- 
lowing this reasoning, one would expect the same 
kind of oscillations also in the AR case for a = 1.0 
(cf. the ordered pattern in Fig 



101. 



In Figs. 15 and 16 the time average polarization index 



The reason 

for not seeing this in Fig. 18 is simple: The average 
distance in e x direction is smaller than the width of 
the support of our induced density. The total density 
therefore smoothens out the periodic structure of the 




FIG. 17. The projection of the pedestrian density on the e x 
axis pg^ as a function of the spatial coordinate r/, relative to 
the center of mass rj crn . Results in the R case for a = 1 at 
time t = 100s. 
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FIG. 18. The projection of the pedestrian density on the e x 
axis p ea . as a function of the spatial coordinate r/, relative to 
the center of mass r] cm . Results in the AR case for several 
different values of a at time t — 100s. 



individuals' positions. What we do see is the fact 
that in this situation the individuals do not occupy the 
whole corridor, but are confined to a certain section of it. 
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FIG. 19. The Morisita index J as a function of the number 
of pedestrians N. Results in the AR case (vertical axis left 
hand side, black) and R case (vertical axis right hand side, 
grey), for several different values of a at time t = 100s. 



(like in the high density core). There is no (sufficient) 
compensation driving the individual forward as this 
effect is decreased/switched off by lowering <j. 
Apparently the mechanism that drives individuals to 
the back has more effect as a decreases. 



In Fig. 18 we 

namely see that due the periodic boundaries the tail 
already reaches again the front of the core group. We 
expect that the Morisita index (see Section III B 3 1 will 
provide extra support for our observations here; this 
index should be high when pedestrians concentrate more 
and more in a particular region. 



3. Morisita index 



What we observe is that including attraction in 
the interaction has some regularizing effect. For each of 
the three choices for a, there is a core of high density 
around the center of mass, which forms the heart of our 
crowd. This core is present due to the initial condition 
that was concentrated on a part of the corridor. Without 
attraction the repulsive interactions would drive it apart. 
However, attraction is not able to completely diminish 
the effect of repulsion. Especially when a decreases, 
the projected density of the core also decreases, while 
there is a tail of mass just behind it (that is, in the 
graph on the left of the center). This effect can be 
explained by the fact that individuals in the anisotropic 
case are driven backwards if they are too close together 



Fig. [19] shows the Morisita index as a function of the 
number of individuals N . The corridor is subdivided in 
M = 6400 boxes. Note that in the graph the results of 
the AR case and R case are combined; each case has its 
own scaling on the vertical axis. 

The value at N = 100 for a = 0.0 in the AR case 
is omitted, since this is very big. This is exactly the 
case of Fig. [8j in which symmetry is no longer present 
in the distribution of the individuals. We expect this 
to be related to the unphysically high Morisita index. 
We should not use the results corresponding to this 
situation for basing our conclusions on. Note however, 
that in Fig. [16] the graph for a = 0.0 does not arouse 
any suspicion at = 100. 
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The increase in the Morisita index as N increases 
(in the R case ) can be explained. For repulsive inter- 
actions, the individuals have a tendency to move as far 
apart as possible (if possible until they are a distance 
Rft apart). This was illustrated by Figs. [HJ^T} As N 
increases, they are however packed together more and 
more, thus leading to an increase in Morisita index. If 
we assume (to obtain an approximate result) that the 
individuals are distributed uniformly, then n,- = j^j for 
all i. Taking into consideration the ordered distribution 
in Figs. [5}[7l this assumption is justifiable. It follows 



M 



M - 1 



that Im = — = 1 — — — > 1 as N — s- oo for 

N - 1 N - 1 

fixed M. Moreover, this implies that Im tends to its 
limit value from below. This matches with the increase 



of the curves in the R case (Fig. 19 1, and we conjecture 



that the Morisita index will tend to 1 for N increasing 
beyond N = 100. 

In the AR case , we observe that the Morisita in- 
dex is not monotonic in a for N smaller than N rs 45 
where the graphs for a = 0.0 and a = 0.5 intersect. 
This point of intersection is hard to explain and requires 
further investigation. Moreover, it is hard to draw any 
conclusion about the precise dependence of the Morisita 
index on N or a in this case. One could argue that 
decreasing a corresponds to lower Morisita index. At 
least this is the case if one compares a = 1 to a < 1. 



IV. CONCLUSIONS AND OUTLOOK 

The behaviour of a crowd of pedestrians inside a corri- 
dor, in which the individuals interact via an anisotropic 
way, can be distinguished clearly from the case in which 
pedestrians interact in a completely isotropic way. In 
particular, we observe the following differences compared 
to the isotropic case: 

1. the polarization index increases with increasing 
anisotropy (i.e. decreasing a); 

2. the projected density along the e x axis shows a 
symmetric profile around the center of mass for the 
isotropic case. However, increasing anisotropy im- 
plies loss of symmetry; 

3. the Morisita index, as a measure of clustering, de- 
pends clearly on the anisotropy. It increases (with 
increasing anisotropy) in the R case and, roughly 
speaking, decreases in the AR case . 

4. In case of repulsive interactions, the crowd tends 
to fill the whole corridor. If attraction is included, 
the group stays compact. Decreasing a however 
seems to diminish this kind of social cohesion as 
individuals do not look behind. 



As a result of the study presented in this framework, 
many new questions arose. Future research should be 
concentrated on the following three directions: 

• Most obviously: what is the effect of a further in- 
crease of the number of pedestrians? Do the ob- 
served relations between the measured quantities 
and the number of pedestrians still hold? How 
do the observed limiting values of the polarization 
for large t depend on A? What about the limit 
N — > oo in this case? 

Can we gain more insight in the strange issues of 
the AR case (cf. Fig. 19 intersecting curves for 



a = 0.0 and a = 0.5; extremely high Morisita index 
at N = 100 for a = 0.0) by further measurements 
of the Morisita index for increasing A? E.g. can 
we extrapolate information for N > 100 back into 
the interval [0,100]? 

If N increases, the natural thing to do is to consider 
the discrete-to-continuum limit (i.e . construct ed- 
ucated procedures to derive mean-field limit equa- 
tions). Does such limit exist, can we derive it, and 
can we compare the effect of anisotropy in the limit 
to the observations of the current work? 

How much does the large time behaviour of the 
crowd depends on the initial conditions? The ini- 
tial distribution of pedestrians in this thesis is not a 
realistic situation. People starting to enter a corri- 
dor are in real life never distributed in a crystalline 
structure manner. However, for an escape situation 
(for example in the case of fire) it seems reasonable 
to assume that a group of people starts, being clus- 
tered, at one side of a corridor. Therefore, as an 
extension of this research, we propose to use as an 
initial distribution a more realistic configuration in 
which people are placed at one side of the corri- 
dor, with their positions slightly perturbed from 
the grid points. Averaging over a large collection 
of such perturbed initial distributions, will lead to 
effective results. Are these averaged results compa- 
rable to the ones presented in this paper? In other 
words: is averaging the results basically the same 
as removing the fluctuations from the initial condi- 
tions? 

With respect to random initial conditions, we are 
also interested in whether these do or do not allow 
order and patterns to be formed (cf. Figs. [5]-{7]and 
10). Do random initial conditions lead to 'preser- 
vation of chaos' ? 



• What happens if we try to make our model more 
realistic: e.g. change the shape of the domain f2, 
or allow variation in the direction and magnitude 
of individuals' desired velocity? Including more so- 
phisticated active parts in the boundary (doors) or 
impermeable objects within the domain, automat- 
ically leads to questions about the efficiency of the 
flow. Which geometry leads to the fastest evacua- 
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tion? First steps in this direction have been made ACKNOWLEDGMENTS 
in 0. 



The issues addressed in this paper show that anisotropy 
related to perception has nontrivial effects on the global 
dynamics of a crowd. Certainly, these effects cannot be 
neglected. More work, both numerically and analytically, 
is needed to extend and formalize our results. 



The authors thank C. Storm, F. Toschi, F. van de Ven, 
H. Wyss (all with TU Endhoven) and R. Fetecau (Si- 
mon Fraser Univ. Canada) for a series of fruitful dis- 
cussions on the dynamics of self-propelled particles with 
anisotropic motion potentials. 
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